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I. INTRODUCTION 

Several laser interferometric gravitational-wave (GW) 
detectors, including LIGO, VIRGO, GEO600, and 
TAMA300 are currently searching for the GWs emitted 
by inspiraling compact binaries, which consist of black 
holes and/or neutron stars. Because the data analysis 
method used by these experiments, namely, the matched 
filtering technique, requires the detector's output signal 
to be compared with a large amount of theoretical wave- 
forms (templates), the post- Newtonian (PN) calculation 
for both of the binary's motion and the gravitational 
waves emitted has to be performed. For nonspinning bi- 
naries, the PN expansion has been successfully carried 
out through 3.5PN order. However, since astrophysi- 
cal observations suggest that most astrophysical objects 
carry a certain amount of spin angular momentum, and 
compact objects like black holes are usually rapidly ro- 
tating, the effect of spin is too large to be ignored. 

Both the dynamics of a spinning binary and the GWs 
emitted by such a system are very different from those 
of a nonspinning system. The coupling between the or- 
bital angular momentum L and the individual spins Sq 
leads to precession of the individual spins and the or- 
bital plane, which in turn leads to additional amplitude 
modulation of the GWs emitted by the system. The 
detailed calculation by Kidder at IPN order [T] showed 
that spin itself can also directly contribute to the gravita- 
tional waveforms and the emission of energy and angular 
momentum. Including spin as an intrinsic parameter of 
the source also increases the dimension of the parame- 
ter space to be used in the data analysis process, which 
not only requires more computational resources but could 
also affect the accuracy on parameter estimation [51 [3] . 
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Thus, it is desired to carry the PN approximation for 
spinning binary systems to a sufficiently high PN order. 

The spin effect to the motion and to the gravitational 
field have been a long-standing problem in general rel- 
ativity (GR). Papapetrou and Corinaldesi in the 1950s 
[H [5] calculated the leading order spin effects to the mo- 
tion of a spinning test body in a given gravitational field. 
The leading order spin-orbit (SO) and spin(l)-spin(2) 
(S1S2) contributions to the equations of motion for a sys- 
tem of two spinning black holes were derived by D'Eath 
[B], Barker and O'Connell [71 H] in the 1970s, and later 
by various authors (for references and reviews, see, e.g., 
[9l|T0]). Using Thome's multipole expansion formalism 
[TP in terms of symmetric trace-free (STF) radiative mul- 
tipolcs in the harmonic gauge, the group of Kidder, Will, 
and Wiseman derived the leading order SO and S1S2 con- 
tribution to the gravitational radiation flux fV, fT2] for a 
general spinning binary system and the polarized gravi- 
tational waveform emitted by a spinning binary system 
with quasicircular orbit The leading order SO and 
S1S2 radiation reaction effects to the equations of motion 
(EOM) were derived by Wang and Will in the harmonic 
gauge [131 [13] . A general form of the SO contributions 
in arbitrary coordinates was also derived by Zeng and 
Will using an energy and angular momentum balance 
approach [15]. 

The (conservative) next-to-leading order (NLO) spin 
effects were only tackled recently. The first derivation of 
the NLO SO EOM was attempted by Tagoshi, Ohashi, 
and Owen |16j . their result was essentially confirmed by 
Faye, Blanchet, and Buonanno in the harmonic gauge 
[17) . and later also by Damour, Jaranowski, and Schafer 
using a Hamiltonian approach in the Arnowitt, Deser, 
and Misner (ADM) gauge [I8 (see also 19]). The corre- 
sponding energy flux (and formulas for a phasing) is given 
in [5D|. The first attempt to compute the NLO S1S2 con- 
tributions to the EOM of a spinning binary system was 
made by Porto and Rothstein [21] using an effective field 
theory technique, namely, an extension of nonrelativis- 
tic general relativity to systems with spin . The 
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first complete NLO S1S2 Hamiltonian was presented by 
SteinhofF, Schafer, and Hergt \W\ [21], and agrees with 

Though the above mentioned results are useful for the 
creation of templates, further work needs to be done. In 
general, a parametrization of the orbits must be obtained 
by solving the EOM. It is common to describe the conser- 
vative dynamics in terms of certain orbital elements, see, 
e.g., [27] • Spin precession and dissipative effects can then 
be described by secular EOM of the orbital elements. For 
explicit solutions including spin see [2H1 [IS]- It is also 
possible to obtain the dissipative orders of these secular 
EOM with the help of the conservative parts as well as the 
energy and angular momentum flux. In this way secular 
EOM corresponding to the LO radiation-reaction EOM 
mentioned above have already been obtained in [5DH33j . 

In this paper, we extend the canonical formalism of 
ADM [34H36] to n-body systems with n spinning objects 
up to a PN order sufficient for the computation of the 
SO and S1S2 contribution to the next-to-next-to-leading 
order (NNLO) conservative Hamiltonian H^^p^^ and the 
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It is im- 



leading order dissipative Hamiltonian H', 
portant to mention that we count PN orders in a rather 
formal way; see Appendix |A] The canonical framework 
in the present paper is constructed order by order in the 
PN setting by utilizing the global Poincare invariance as 
the important consistency condition, similar to [19| . Fur- 
ther a general formula for the interaction Hamiltonian 
between matter and the transverse-traceless part of the 
metric hj"^ at 3.5PN is derived. From this Hamiltonian 
a wave equation for hj^"^ can be obtained by canonical 
methods. For the canonical formalism presented in this 
paper, this wave equation agrees with a corresponding 
one that can be followed directly from the Einstein equa- 
tions. This provides a thorough check of the canonical 
formalism. Further the obtained formulas are the basis 
for applications. Using the wave equation, we are able to 
derive, in the radiation zone of a system with two spin- 
ning objects, the leading order SO contribution to hj^^ 
and the energy flux (see also U [33 ) • 

A canonical framework for spinning test-particles valid 
to any PN order and linear in the spin of the particle was 
given very recently in |38j . The Hamiltonian of a spinning 
test-particle in Kerr spacetime was given explicitly. This 
includes parts of the conservative Hamiltonian -ff<°3PAr 
mentioned above (as well as contributions of cubic and 
higher order in spin; see also [33 HD])- An action ap- 
proach to the canonical formulation of self-gravitating 
spinning objects valid to all orders linear in the single 
spin variables was recently given in 41 and is shown to 
agree with the present paper up to 3.5PN. The order by 
order construction performed in this paper gives an in- 
dependent derivation of the results in [IT] up to 3.5PN. 
Further the method developed in this paper to construct 
a canonical formalism might have some advantages over 
an action approach at higher orders in spin. Knowledge 
of the formalism in [19] and its extension given in the 
present paper was important information to succeed with 



the action approach in [JT]. The consistency checks and 
applications performed in this paper also apply to the ac- 
tion approach. Further, the present paper provides more 
details than [IT]. In particular, the source terms of the 
constraints depending on canonical variables are given 
explicitly. 

The paper is organized as follows. In Sec.|TT]it is shown 
how the ADM formalism can be extended to spinning ob- 
jects order by order in a PN setting. In Sec. [Hi] this ap- 
proach is applied to 3.5PN and linear in the spin. A com- 
parison with the action approach in |41j is given. Section 
IV gives the PN expansion of the constraints, including 
the matter source terms in canonical variables. In Sec. 
[V] a general formula for the interaction Hamiltonian is 
derived. Further it is shown that the evolution equations 
given by this Hamiltonian, specialized to our canonical 
formalism, coincide with the ones following from the Ein- 
stein equations. In Sec. |VI| the SO energy flux is com- 



puted. Section VII gives conclusions and outlook. 

Our units are c — 1 and G — 1, where G is the Newto- 
nian gravitational constant. Greek indices will run over 
0, 1, 2, 3, Latin indices from the middle of the alphabet 
over 1,2,3. Latin indices from the beginning of the al- 
phabet label the individual objects. For the signature 
of spacetime we choose +2. The short-cut notation ab 
{— ai^'b^ — a^b^^) for the scalar product of two vectors 
and 6^ will be used. Square brackets denote index an- 
tisymmetrization and round brackets index symmetriza- 
tion, i.e., a^'^fe'') = i(a''6'' + a'^b''). The spatial part of 
a 4- vector a; is x. Round brackets around an index de- 
note a local basis, while round brackets around a number 



denote the formal order in 



as m 



II. FROM ADM ENERGY TO ADM 
HAMILTONIAN 

In this section we outline how the ADM canonical for- 
malism [34H36j can be extended to self-gravitating spin- 
ning objects order by order in a PN setting. We only con- 
sider a fully reduced canonical framework here where the 
gauge is fixed and all constraints are eliminated. Then 
the ADM energy can be used as a Hamiltonian, the ADM 
Hamiltonian, if it is expressed in terms of variables with 
standard canonical meaning. The transformation to such 
variables can be found from consistency considerations. 
At the 3.5PN SO and S1S2 orders the global Poincare 
algebra and the constant Euclidean length of the canon- 
ical spin uniquely fixes this transformation to standard 
canonical variables. 

The approach outlined here is a natural generaliza- 
tion of the one in [19]. It was suggested in Appendix 
B of [19 , by considering the algebra of the gravitational 
constraints, that at higher orders spin corrections to the 
canonical field momentum might be necessary, and that 
the gauge structure needs to be extended. Indeed, the 
former is an important ingredient of the approach in this 
paper [see Eq. (2.6|] as well as of the action approach 
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in [IT]. Further, the action approach is based on tetrad 
gravity, which has more gauge freedom than metric grav- 
ity. In this paper, however, the gauge is always fixed and 
the original gauge structure is less important. 



Field constraints 



by 



' can * can,_7 



Vi 



"ij *^can,fe 7 



Tt' ■ +TT-' 
can, J ' "can.z 



3 

r, "tj "can,fc 



2 ca.n,ijk ' 



(2.10) 

jk ■ 
(2.11) 



Most important for an explicit calculation of the gener- 
ators of the global Poincare algebra, including the Hamil- 
tonian, are the constraint equations of the gravitational 
field. They can be written as 



1 



7R^ 



1 



J- jk 

-8^7.,- .,k 



lijlklT^^'^T^^^ 



(2.1) 
(2.2) 



It holds 



"can ^ can J canj 



_ijTT _ cTTij kl 
''can ~ "fc; ''^can i 



(2.12) 

(2.13) 
(2.14) 



with the inverse Laplacian A ^, the partial space- 
coordinate derivatives di and 



with the definitions 



rTTfc/ 



matter 



-V7(7^'7^' ~ 7*^7'')^fc/ , 



(2.3) 
(2.4) 
(2.5) 



and arise as certain projections of the Einstein equations 
with respect to a timelike unit 4-vector with com- 
ponents = (-iV, 0,0,0) or {1,-N')/N. Here 
"fij is the induced three-dimensional metric of the hy- 
persurfaces orthogonal to n^, 7 its determinant, R the 



--l[{6u-A-'dA)i6jk-^~'d,dk) 
+ {6,k-A-^d,dk){Sji-A-^djdi) 



(2.15) 



See Sec. IV for details on the decompositions for 7^.,- a nd 



T*-^„. Notice that the form of the trace term in (2.8) is 



adapted to the Schwarzschild metric in isotropic coordi- 
nates, with obvious advantages for perturbative expan- 
sions. 

Now the four field constraints can be solved for the four 



three-dimensional Ricci scalar, K,, the extrinsic curva- variables and 7r*^„ in terms of hj^ , ttJ^J"^ and matter 



ture, N the lapse function, TV' the shift vector, y^T^^^ 
the stress-energy tensor density of the matter system, 
and ; denotes the three-dimensional covariant derivative. 
Partial derivatives are indicated by a comma. 

For nonspinning objects 3^ is the canonical mo- 
mentum conjugate to jij before gauge fixing. For spin- 
ning objects we now make an ansatz for the canonical 
field momentum of the form 

■""can = TT -|- TTj^atter i (2-6) 

where Tr^atter shall be linear in the spins and will be fixed 
later on. In the ADM transverse-traceless (ADMTT) 
gauge defined by 



variables, which enter through the source terms T^™^,**^'' 
and 'H™'***"''. An analytic solution for (j> and TTcani how- 
ever, can in general only be given in some approximation 
scheme. 



In the ADMTT gauge the momentum constraint (2.2) 
can exactly be written as 



- «7 

TT ■ — 
can. J 



can, J 



tcr _|_ -^Trn 



(2.16) 



l_jfcTT, TT _ (---jfcTT. TTn . 
r)"can jk.i V'can "fci J ,J J 



with the definitions 



TT 



(2.7a) 
(2.7b) 



which will be used throughout this paper, one has the 
decompositions 



TT + TT 

" ran ' " car 



(2.8) 
(2.9) 

TT _ 



where hj^^ and 7r*{'J'^ are transverse-traceless, e.g, fi. 
hj^^j = 0, and 7rj{„ is related to vector potentials T4an ^^'^ 



can 



1-(1 (Tf, 



can 



^ "can J 



+ KanC'-feij' + f^kl,i ~ f^ij^k) 



TT 



r TT > 



2 can,k ij ' 



(2.17) 

— [-2(7,fc7r^4tt„)j -I- 7r^attcr7j7c,j] ■ (2.18) 



This equation will allow us to derive explicit expressions 
for total linear and angular momentum without solving 



the constraints. Notice that 



Ail, and All^ = 0. 
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B. Global Poincare algebra 

The global Poincare algebra is a consequence of the 
asymptotic flatness and is represented by Poisson brack- 
ets of the corresponding conserved quantities. These 
quantities are the ADM energy E, total linear momen- 
tum Pi, total angular momentum Ji = \eijkJjk, and the 
boost vector . They are given by surface integrals at 
spatial infinity. The boosts have an explicit dependence 
on the time t and can be decomposed as = — tPi, 
where X* = /E is the coordinate of the center-of-mass. 

will be called center-of-mass vector in the following. 
The corresponding surface integrals read, with spatial co- 
ordinates denoted x*. 



E 

Pz = 
J a — 



1 

1 

16^ 
1 
1 

'8^ 



(2.19) 

(2.20) 
(2.21) 

(2.22) 



See, e.g., Using the gauge con dition s and also the 

momentum constraint in the form (2.16), these surface 



integrals can be transformed into the volume integrals 



E = 


WttJ 


^ d^x A(j) , 


(2.23) 


G = 




^ d^xx'Acj), 


(2.24) 


Pz - 


pmattcr 


16n / can "■fci,i J 


(2.25) 


Jij = 


Tinattcr 
'Jij 


-±-jd'x2{.fI'hlJ 


jfeTT, TT^ 
"can "•fei ) 



1 

167r 



{-"^ ~can "-fc(,7 ^ "can ""fc/.J : 



with the matter parts 
P/"^""' = J d^x {n^ 



Jf. 



/ matter 



/TTinatter 



x^m 



x'n 



^j^^mattcr) 



(2.26) 

(2.27) 
(2.28) 



Here we used the fact that has a compact support. 

Now we require that the matter parts of total linear 
and angular momentum are of the form 

^matter ^^p^^^ (2.29) 

a 

jrnattor ^ ^(^^p^^. _ ^Jp^^) + ^ gai^)^,) , (2.30) 



where z*, Paj, and 5'a(i)(j) = ^ijkSa{k) are the canonical 
position, momentum, and spin of the particles, because 
this is the expected form for standard canonical variables 
with equal-time Poisson brackets 



{zl,Pa,} = 5.. 



} ^ijkSa 



{k) 



(2.31) 
(2.32) 
(2.33) 



zero otherwise, where e^fe = ^{i — — k){k — i) . Equa- 
) and ( 2.30[ ) ensure that a great part of the 
gebra is fulfilled; see Appendix IS] With the 



tions ( |2.29 
Poincare a 
definition 



Matter = IStT^TT;; 



(2.34) 



where Sa — <5(x — z^) with normalization J d'^xSa = 1, 
the source of the momentum constraint 7^™'^**°'' then is 
of the form^ 



S, 



E 

a 



1 



^T^a iT-kj - iT-kt ■ 



, (2.35) 

(2.36) 



At linear order in spin and 3.5PN, the first equation de- 
fines the canonical momentum, while the second one fixes 
the canonical position, the triad (i.e., the local basis of 
the canonical spin) and tt*-' , up to canonical transforma- 
tion. Further the Euclidean spin length Sa given by 



2si 



'2Sa{i)Sa(i) 



Sa{i){])Sa{i){j) 



(2.37) 



has vanishing Poisson bracket with all quantities, includ- 
ing the Hamiltonian. Thus Sa must be a constant of 
motion. 

The ADM Hamiltonian Haum results as 



^^ADM = / d xA(j)[Z^,Pai,Sa(i),hij , TT, 



can J 



(2.38) 

This is the ADM energy depending on the canonical vari- 
ables. It arises from solving the constraints for 0, once 
the source terms of the constraints, y,™^*-*-'^' and 7^™attcr^ 
are expressed in terms of the canonical variables. An ac- 
tion corresponding to i/ADM is given by (4.33) in [TO] or 
(51) in j41J. 

Total linear and angular momentum could, of course, 
also be represented on the phase space in a more com- 



plicated way than given by (2.25), (2.26), (2.291, and 



(2.30). However, the ADMTT gauge manifestly respects 
the Euclidean group in its standard representation, which 
implies that its generators Pi and Jij are also in its stan- 
dard representation on the phase space; see Appendix [B] 
and also 1421. 



^ Here we assumed that the variables from difFeront objects do not 
mix (e.g., as in Pi<52) at this stage. 
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III. THE SOURCE 



A. (3+l)-split 



The stress-energy tensor density to linear order in spin 
is given by 



dr 



(3.1) 



in the covariant spin supplementary condition (SSC) 

5r^ia. = 0. (3.2) 

Here TOq is the mass, Ua the 4-velocity, r the proper 
time parameter, Sf^"^ the spin tensor, || denotes the four- 
dimensional covariant derivative, and (5(4)a = S{x — Za) 
with normalization J d'^x5(i)a = 1- the coordinate 

of the a-th object. The matter EOM, i.e., the Mathisson- 
Papapetrou equations |31I351|35], in covariant SSC and 
at linear order in spin can be followed from = as 



DS^ 
Dt 
Dp 



0, 



^ 9 pPa"-a'->a J 



(3.3a) 
(3.3b) 



"^g — ,M _ fa 

dr nia 



(3.3c) 



Here '^^^R-^'pAQ four-dimensional Riemann tensor 

and D the four-dimensional covariant parameter deriva- 
tive. The spin length Sa is given by 2s^ =^^Sa^i.v and 



obviously is a constant of motion due to (3.3a). 

The (3+l)-split of ul = -1, the SSC and the spin 
length reads 



2sl 



npa 



(3.4) 

(3.5) 
(3.6) 



Notice that nS^ = 'j^^nSaj- The components of the 
stress-energy tensor density are given by, with Sa = 
'5(x-Za), 



n 



matter ^ ^ 



npaSa - K"'— 5a " {nS'^5a);k 



npa 



(3.7) 



itter _ 

a 

-7^. = E 



c{k nmPaknSaJ 

npa 



PaiPaj , j^k , Pa{iSaj)kPalK''^ nSakPa{tKj^ PaknSa(tK^^ Pa(inSaj)PakPalK 



kl 



+ J 



I Sak{iKj\ I , ,„ 



("Pa)' 



(3.8) 



npa 



-6a - r 



{npaY 



Sa 



where Tij — ^Tij. After transition to Newton- Wigner (NW) variables 

Pai'^Saj I PajlT-Sai „ Pak'J^'^ Saji 

noa 



Saij — Saij 



a a 



ma - npa nia - npa 



ma - npa 

Pat = Pat - nSaKik - TT^'^Tjfc.i 



I m _ PamPg,!^' M.r^k) PapPa^^"" 
o ' li 7D /™ ^ T> \ ' li. /™ _ : 



nPa{ma - nPa) 



ma{ma - nPa) 



(3.9) 

(3.10a) 
(3.10b) 
Sajk , (3.10c) 



where F^^- are the three-dimensional Christoffel symbols, the source expressions of the constraints read [now Sa 
S{x - Za)] 



E 



nPaSa 



Pal 



ifSauPc 



a] I „,mn 



2V nPa 



+ 1" 



S ami Pa j -Pan Pal ^Pal'^aij Pai^aj c- \ r 

{nPaY{ma-nPa) nPa nPa ' 



iria - nPc 



rY'l^'SajkSa - nPaSzlSa^ 



(3.11) 
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E 



P 



PalPak 



nPa{ma - nPa 



(3.12) 



Obviously (3.12) is now of the from (2.35), which 



uniquely fixed the relation between covariant linear mo- 
mentum pai and canonical momentum Pai, Eq. (3.10c I 



(3.6) into 



The spin redefinition (3.10a) transforms the spin length 



250 = "/'''"f'^ SaklSaij = Sa{i)(j)Sa(i)(j) 



(3.13) 



The symmetric part is thus fixed. The antisymmetric 
part e*(^\ however, must be imposed, as it represents the 
three rotational degrees of freedom left in the definition 
of the local basis. Thus e*^-'-' represents the degrees of 
freedom left in the definition of the canonical spin vari- 
able. 



where Sa{i)(j) are the components o f Sgjj in some local 
Euclidean basis. Comparing with (2.37) suggests that 



is equal to the canonical spin Sa 



WO)- 



However, 



a local Euclidean basis is only unique up to a rotation. 



Fortunately Eq. (2.36) will be seen to uniquely fix a basis 
such that and Sa,(i)(j) can be identified in that 

basis. (Notice that Sa{i)(j) are not the components of 
the covariant spin Saij in a local basis here, in contrast to 



[41] .) The redefinition of the position ( 3. 10b ) consists of a 
term known from flat-space and a ye t unk nown quantity 
dzl, which will also be fixed by Eq. (2.36) later on. 



B. Triad 



C. Fixation of the NW variables 

Whereas the canonical mor nentum was already unam- 
biguously fixed by ( |2.35[ ) as ( |3.10c| ), Szi, e''-J\ and <J 
are still unknown. We will see now that these can be 
fixed up to a canonical transformation with the help of 



(2.36). For our source one gets for the leading orders of 



*j - 

"(3) - 

AiJ] _ 
5a(5) - 



S. 



(2)' t;-(2^ Oaik){i) 



PnJ: 



a{2) 



(3.18) 
(3.19) 



The relation between Saij and Sa(i)(j) can be written 
with the help of a triad e*'-'-' as 



(3.14) 



This spin has a constant Euclidean length for all choices 
of e*^*^*^ . Notice that the triad is needed only on the world- 
lines and not as a field over the entire spacetime here. 
The triad can be split into symmetric e*^^^ and antisym- 
metric e''^-'^ parts as e'*^-'^ = e^^^^ + e^^J\ Perturbative 
expansion of e'^'^^e-''-'^-' — 7*^ leads to^ 



and thus from (|2.36|) one concludes 
0, 



tik) 



-(2) 



Sz. 



a(2) 



(3.20) 



It is crucial that e^^^^ must be antisymmetric. At the 
next order it holds 



^a{7) - ^(4) '^a(k){j) (^(4) Ja{k){i) 



(3.21) 



(n—m) (n) ' 



(3.15) 



where 7^^^ = Sij and e^g^'* = Sik was assumed. For exam- 



ple, the leading order results are: 

6(4) 



6(4) - 96(2) 6(2) - J5^3<P{i) + ^^^■'^(2) 



(3.16) 

1. 



, TT 



(3.17) 



^ Notice that in this formula the subscripts in round brackets de- 
note the formal order in c~^, not an index in a local basis. 



Further T^^a(3) ~ ^ because tt*^ is linear in spin, symmet- 
ric, and contains an even number of momentum variables 



(because of parity). Equation (2.36) then leads to 
e;i?=0, Sz. 

M 



a(4) 



0. 



For sjj'^g^ one has 



a(9) 



.i{k) 
-(6) 



Sa{k)(j) 



-(6) 



Sa(k){i) 



Pm5: 



'a{6) 



-Paj'^^a(6) 



TT 



2m2 



2rn? 



PakSa{k)(l)Pa[ihj]l 
PakPalSak[ihJ-^ ■ 



(3.22) 



(3.23) 
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The most general (sensible) solution of (2.36) at this or- 
der is 



jj ^ 1 — C 



{PaiPakSa{k){j) + PajPakS a{k)(i)) , 



'(6) 



c 



2^2 PakPa[ih]yi , 
C 



(3.24a) 
(3.24b) 



^^l{6) ^ ^^PcLjiSa(k)(t)hjk + Saik)U)h,k ) , (3.24c) 

with an arbitrary constant C. Notice that tt^^-j^ = 0. 
Now we can remove the ambiguity C by a canonical 
transformation generated by 



C 
Arnl 



PaiPakS a(k 



mi) 



(3.25) 



which transforms an arbitrary phase space function A as 

A~^A + {A,g}, (3.26) 

to the required order. For our fundamental variables this 
means 



TT 



(3.27) 



AttC 



kl rn'2 ^a-kPamSa(m)(l)5a , 



(3.28) 



6(6) 



'aik){j) - Sa(i)(k) , (3.29) 

(3.30) 
(3.31) 



~^ 

Pn.i Pfi' 



Sz. 



a(6) 

c 



TT 



PalPa]Sa(j)(k)hkl, 



This indeed removes all terms depending on C from the 
source expressions in (3.11) and (3.12) at the considered 



order. We can therefore choose C = 0. This choice has 
the nice properties that the triad fulfills the gauge con- 
dition e**^^-' = e-'^'-' and that the transition to the NW 
position can be expressed in terms of a Lie-shift; see 
Appendix [C] Further this choice agrees with the action 
approach in [¥T] . 



D. Comparison with the action approach 

For the action approach in [41] the gauge condition 
^i(i) ~ ^3{i) ~ ^ij holds to all orders, which fixes Cij as 
the matrix square-root of the three-dimensional metric 
e-ije^k = lik, or 



Vilki) ■ 



(3.32) 



Furth er, th e transition to NW variables in [5T| agrees 
with (3.10) for Sz^ = to all orders. Notice that Eq. 



(31) in [IT] is written in terms of the covariant spin and 



one has to use Eq. (45) in [41] and (3.33) on the triad 
terms. It is shown in [51] that a suitable choice of ttJ* 
extends the ADM formalism for spinning objects to all 
PN orders linear in spin. According to ^Ij, it holds 



1 



where the quantity B^- is defined by 



2 
1 

kl 



^kl-^a 

maPaanS, 



nPainia - nPa) 



(3.33) 
(3.34) 

(3.35) 

which must be evaluated perturbatively using Eq. (4.22) 



^3 

^k[i^j]k,^ ' 

This can also be written as 
25*^' - e -^^ 



Bf.lkL^ 



den 
dlk 



e.g.. 



2B^l = Us.^khU - 5.,,khl]) + Oiih^^r) . (3.37) 



Further it holds 



TTlaPainS, 



2nPa{ma - nPa 



(3.38) 



which follows from B^j6ki = 0. This leads to a deviation 
of our gauge condition 7r*!^,j = from the original ADM 
one tt" = at the formal 5PN order. 



We will now show that (2.36) exactly holds for the ac 



tion approach. This gives a check of the action approach 
to all orders, for any approximation scheme. From the 
discussion above it is clear that the source of the mo- 



mentum constraint is of the fo rm (|2.35|). Further it holds 

si' 



Jik^a'' ■ Equation (2.36) then leads to the condition 







5a(fc)(/) 



2e 



(cfeie 



CkjC 



li 



'2SkiSt 



Evaluating 



results in 
der 



dip 



dip 



-pi<^qj 



(3.39) 



(3.40) 



(3.41) 



The second term on the right-hand side is symmetric in i 
and j and cancels from (3.391 with (3.36) inserted. It is 



then easy to see that (3.391 and thus (2.361 are fulfilled 



The action approach shows that (3.11) and (3.12 ) can 



be applied to all orders, if Szl — 0, gHlT— q, and (3.33) 
are inserted. Further, if the spin correction to the field 



momentum tt^J is neglected, (3.11) and (3.12) coincide 
with (4.23) and (4.25) in [19]. In |19| a Lie-shift was used 
to redefine the position instead of a Taylor expansion 
here; see Appendix [C| 
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IV. PN EXPANSION 

In this section we give the PN expansion of the field 
constraints relevant for the ADM Hamiltonian up to and 
including 3.5PN. Here and in Sec. |VI| we made use of 
xTensor [17] (a free package for Mathematica [IH]), espe- 
cially of its fast index canonicalizer based on the package 
xPerm [15] . 

For TT*^ we use the general decomposition 



with 



Sir'-'' . 



kk 



(4.1) 

(4.2a) 
(4.2b) 
(4.2c) 
(4.2d) 



Notice that it holds nH^ ^ = n\j + 7rli^^^^^. j and thus 

"can ^ " ^'^j'-^ "matter- 

The expansion of the momentum constraint immedi- 
ately follows from the exact formula 



1 



i/cTTtTT / ifcTT7TT\ 
2^ 'T'jk,i - (t^ hz ),j : 

1 



which is analogous to (2.16). With the help of 



(4.6) 

(4.7) 
(4.8) 



This can be shown by inserting (4.2 ) and (2.151 into (4.1 ), 



which then turns into an identity. We further introduce 
an alternative vector potential by 



for which it holds 



1 



(4.3) 



(4.9) 



the expanded momentum constraint can be solved itera- 
tively for tt* by applying an in verse Lap lacian to it. tt*^ 
and then follow from (4.2b) and (4.3). The expansion 
of the source -H^attor given by ([2.35| and 



(4.4) 



Notice that this decomposition reduces to the one for 
7 r^{„ , as 7r^{„ — follows from the gauge condition 



2.7b) 



4.2d 



This gauge condition can also be inserted into 
in the form tt" = — tt^j^jj^j.. This immediately 
yields tt*-' , which is zero at the considered order due to 



^(5)mattor ~ (3.38) for higher orders. Similarly 

the decomposition of can be derived using (2.7a I. 

Now we have to decide whether to use i^"^^^ and 7r^{jj, 
or TT*^"""^, tt'-' , and tt*-' for the expansion of the field con- 
straints. We choose the latter option as it simplifies the 
calculation at the considered PN order. Then one has to 
go over to 7r*{J'^ later on using 



" "can 



"fc; " matter 



(4.5) 



*a(3) ~ '^a(i)O') ' 



S'!,^-. = -^^PakPaSa{]){k) + (« j) , 



2to^ 



(4.10) 
(4.11) 



3P2 

o'J - a P , P ^ 

A ak-^ ai^a 



Arri^, 



\hHSa(j){k) + i) ■ 



PakPaiSa{j)(k)4>(2) 



(4.12) 



Notice that 



■a(9) 



is not needed for the Hamiltonian at the 



considered order, it only contributes to total linear and 
angular momentum. 



The expansion of the Hamilton constraint (2.1) reads 



1671- 

1 

16^5^ 
1 

16^ 



A0(4) 



-T/ matter 
rt(2) 

-27 matter 



" 1/ matter i 

:«(2) 0(2), 



-T/matter i2 
(3^ ^(2) <^(2) 



1 

1(3^ 



(4.13) 
(4.14) 

(4.15) 
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1 

1677 



-T/mattcr 

- "(8) 



IGtt 



matter 
(4) 



/matter 
t(2) 



512^(2) 



matter j 3 



1 



2^(3)^ 



2) ^ 



1 

16^ 
1 



70(2) (tt'-' 



(0(4) /l. 



ij ) 



0(4) + ^ 

2 .... 
(3)y +^'^(3)^5) 

1 



0(6) 



+ ^ (^(4) 



matter 

^(2) 



1 

16 



■■+ A 

,v 4 



(0(2) 0(2) J /ly'^) 



-^A(/.-^^ 



(2) 

1 / , TT \ 2 

4 
1 



TTi^TT\ 



2 ('^ij "ifc y jfc 



A 1 o /matter /matter i , ^/matter i . ^/matter i , o /matt 

A0(io) = rt(io) - g (^rt(8) 0(2) + rt(6) 'f'(4) + rt(4) 0(6) + ^(2) 

+ ^ (^r6)"">(2) + 2H'('lf *^^0(2)0(4) + 2Hf2f *"-0(2)0(6) + H^^f '"0^4)) 



+ 



64 
1 

512 
1 

16^ 



^(ir'''^(2)+3^(2)""'4)0(4) 



'T/ matter ^4 

4096 (2) 



.4 T/mattor /'/jTTN 



^ (0(4) (^3)) +20(2)^('3)^|'5)) + ((^(5)) 



^'^C3)'^(7) 



(2)^^)-^^"" 



+ ( -^0(4),»0(2),j + 4^0(2)0(2),«0(2)j + ] hj^ 



7 

32 



0(2) e^U-^fe) 



128 
1 

16 



0(2) (/ir/l^^ 



(td) 



(4.16) 



(4.17) 



where (td) denotes a total divergence. These equations can be solved iteratively for by applying an inverse Laplacian 
to them. The ADM Hamiltonian (2.381 results from an integration over the right-hand sides of these equations. The 
source expressions are given by 



n 



(2) 



'ly matter 
rt(4) 



o/ matter 
«(6) 



'T/mattcr 
«(8) 



^rUaSa , 

a 

E 



p2 I 



(4.18) 
(4.19) 



E 

a 

E 



[ (P^)% 


p2 

a 


0(2)^a - 


1 

Arua 


X'- a/ ^ 1 


(P^)' 


0(2)^a ^ 


5P^ 


16m^ " 


8m3 


64ma 



8m3' 



4ma 



2m„ 



(4.20) 



1 



1 



f'«S'a(i)(j)0(2)0(2),j<5a + Tf--P«5'„(,)(j)0(4) j(5a + 7^ PatS aU^hJ^k^a 



E^. 



32TOa 
(P^)^ 



p2 5 2. 

lg^5"^"'^'^(i)(j)'^° + ■^P^^^a{€)U)(l)(2)5a + -^^ParSa{i)U)(l}f2)^a - -^ParSa(i){j)(l>(i)Sa 



-PaiSa{k){i)hJiF 6a Pn.iS„ih\( i\hilJ' 5„ 



4mn 



ai^a{k)(j)'Hk "a 



(4.21) 
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'17 matter 
'^(10) 



E 



5iPir 3(p^)3 

TT^ ^<P(2)0a 



128ml 



1 



256TOa 
r.TT 



32mQ 



3(P^)^ 



(I ^ (2 



9p2 

^f'«S'a(i)(i)'/'(2)0(2)j'5a 



15 



64m^ 
5 

32ma 
1 



256TOa 



^'MS'a(i)(i)0(2)0(2)j<5a - ^^-P«5'a(^)(j)<?!'(4) j<5a 



1 P2 

PazSa{i)(j) ((/)(2)0(4)) j<5a + ^-P«5'aWO)0(6) jt^a - ■^PatSa{j)(k)hJ^k^a 



3 1 

2^^P<^i^aU){k)(l}{2)hJ^kSa + ■^^PaiSa(j)(k)(l)(2)^jhJ^ 5a - -^^Pa^Sa(^){k)(t){2),jhJ^ Sa 



(td). 



(4.22) 



Although the source terms given in this section seem 
to be SO couphngs only, the expressions given here are 
enough to give all S1S2 contributions, too. All S1S2 terms 
in the Hamiltonian come in from the nonlinearities on 
the right-hand sides of the expanded constraints, as in 
[ini m]. At higher orders in the single spin variables, 
however, more contributions are needed in the source of 
the constraints (i.e., in the stress-energy tensor). 



where hJ^'^'^'^^^ denotes the parts of the ADM Hamilto- 
nian depending on hj^"^ and Tr*^^'^. The field EOM then 
read 



Idtt -' 



167r 



Sir' 



d SH 



int 



1 

16^ 



(5.4a) 
(5.4b) 



FIELD EVOLUTION 



TT 



In this section we derive the wave equation for h, 
from the ADM Hamiltonian and compare it with the 
corresponding one that follows directly from the Ein- 
stein equations. Agreement is found, which proves that 
the constructed ADM Hamiltonian gives the correct time 
evolution for the gravitational field up to and including 
3.5PN. This provides a thorough check of the canonical 
formalism derived in this paper and also of Ref . |41| . 



with O = A — and the partial time derivative dt- 

Notice that 0(6), 0(8)j and Tr^:^^ depend on hj^^ and/or 
TT*-'"'""'", and at a first look it seems that one has to ex- 
plicitly solve the constraints for these functions in order 
to get the interaction Hamiltonian. However, one can 
use the expanded constraints to eliminate 0(g) , 0(8) , and 



by performing certain partial integrations. With the 



definitions 



A. Interaction Hamiltonian and wave equation 

The field EOM can be obtained from the ADM Hamil- 
tonian by 



'1(4) 



»2(4) 



— 1 17 matter 



-16TrA" 



(4) 



1/ matter i 



(5.5) 
(5.6) 



1 ; TT _ x'^'T^'-j SHabm 



1 



•jjTT 



can 



167r 



(5.1) 
(5.2) 



where the dot denotes a partial time derivative. It is 
suitable to introduce an interaction Hamiltonian H™^ be- 
tween matter and gravitational field as 



rrint TrTT-parts 

— -"ADM 



1 

16^ 



{^j^kY + ("'can'^)^ 



(5.3) 



for which 0(4) 
grations read 



1/ matter ^ 

"•(2) <; 



1/ matter 
«(4) 



(6) 



'H(^)"'"'0(2)0(6) 



'1(4) 



1 

32^ 
1 

32^ 
1 

47r 



2(4) holds, these partial inte- 



[0(2),»0(2) j/^i^-^] + • • ■ , (5.7) 

[01 (4), » 0(2) J ^i^-'^] H , (5.8) 

[02(4),i0(2)j/li^-'^] H , (5.9) 



11 



n 



(2) 9(8) 



(8) 



l^jmattcr /'^,TT^2 , 1 
167r 



2) + 9«(2) i'^^J j + 



TT 



16 



(5.10) 



_ 1 p,~\li -T/matter 



where dots denote total divergences and/or terms independent of hjl^ and tt]?^ . The interaction Hamiltonian then 



results as 



H 



int 



■-(2) 



matter/ can 



(5.12) 



where 



matter 
(8) 



Sh 



TT 



(6)y 



167r 

5 / \ ■ ■ 1 ■ 

+ ^0(2)0(2),»0(2) J + 2^^^3) [^fg)^, - ^^3)^fe j + '^^"(iykKi) + 2^(3)''f3).fc ' 



The field EOM finally result from ([5^ as 

2B 



□ /jTT _ ^TTij 
'-'"'ij — "fci 

1 • 



(4)fe/ + '^Bf^Q^i^i 



16n 



TT 



dt 



V(3)0(2),; 



A;/ 

(5)mattcr j ' 



with B(6)ij i3(6)ij + TT^^jjj^^^^^^. For our source, one gets 



B, 



B, 



16^ E 

a 

16^ E 



2m ^'^^^'^i^'^ 2rn ^'^'^^<^U)W^°'''^ 



(2),^'?'(2)J : 



p2 5 p2 



4m3 



4m3 



PalPajPakSa(l){i)Sa,k 



5 11 

+ 8^^"'^'^(J')('=) ('^(2)^'»).fc + ^^«'5'a(fc)(i)<?^(2),/c'^a - ^^PakSa(km(l){2)Ja 



+ \Sa(k)(i) (V(^3)^fe + V(^3) j) 



1 3 5 ■ / 

+ ;,01(4)<?!>(2),»j + g02(4)0(2),zj + ^<?!'(2)0(2),»0(2) J + 2^f^3) [^i-|=3)_, " 



64' 



(3),fc 



1 



(5.13) 



(5.14) 



(5.15) 
(5.16) 

(5.17) 



^z-ij \rk I ^ ~ij ~k 
'^^{Z),k^{'i) + 2 (3) (3),fc ■ 



(5.18) 



r 



Here the 0(g) terms in 'H^Jg**'"' were rewritten as quite general and in principle applicable not only to lin- 
_i^matter^^g^ + (td) and handled by jsisl. The time ear order^ spin. Further, for nonspinning objects the 



derivative ''■(s^jj^^ttcr calculated using leading-order 
Hamiltonians. Notice that the formula given for the in- 
teraction Hamiltonian (and thus the wave equation) is 



result in [50] is reproduced. 



One can remove 7r!is , . from the interaction Hamil- 

(ojmattcr 
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tonian by a canonical transformation generated by 

[ <^^<m^torh^^ (5-19) 



corresponding to the choice C = 1 in Eqs. (3.241 



while 7ij^ 



leads to 



with the solution 



(5.26) 



B. Comparison with the Einstein equations 

The time evolution parts of the Einstein equations in 
the variables used here read 

= 2iV7-i/2(^^^. _ i7,,7M7r'=') + N,., + N,,, , (5.20) 



2 'Jmn 



.rnn tj 



i7V7-i/^^^(7r""^™„ 
2iV7-i/2(^^^^^^-V"-'' 



(5.21) 



(Notice that there are misprints in Eq. (8.4) in [in]-) Af- 
ter using constraints and coordinate conditions, this can 
be compared to the results of the last section. Lapse and 
shift are fixed by the requirement that the gauge con- 
ditions (2.71 are preserved in time. From J^jTr*^ = it 
follows 



A7V(2) = In 



(2) 



AiV(4) -4^7^4)ri+4^H^f°'■ 



7^•H 



(2) 



1 

32 



(5.22) 
(5.23) 



with solutions 



N, 



(0) 



N, 



(2) 



(4) 



47rA-i7; 



(5.24) 

1(4) - ^02(4) + ^0(2) : (5-25) 



N, 



(3)i 



-2V, 



(3) 



(5.27) 



Using the Hamilton constraint and these expressions 
for lapse and shift, the PN expansion of the TT-projected 
evolution equations reads 



/^r=2-™^'+2c^^'(^(3)^(2).), 



-Ah. 



inn 



TT 



"kl 



matter 7 TT 
(2) "-kl 



(6)kl 



TT 

kl.m 



with 



(5.28) 
(5.29) 

(5.30) 



B{6)ij — -87r7(6)ij + 107r7(4)ij(?!'(2) - 27r0(2)_y A ^T(4,)kk 
1, ^ 3, , 

+ T01(4)0(2),ij + g<P2(4)<P(2),jj 



64 



0(2)9!'(2),,?!'(2) J + 27f^3) (^fs), 



'(3),fc 



■ 27r(3),feV(3) -f -7r(3)7r(3)_fc . 



(5.31) 



This can be written as a wave equation identical to ( |5.15[ ) , 
but now with different expressions for -B(4)ij and i?(6)ij. 
Also, Eq. (5.28) obviously is the same as (5.16). The 



question thus is if the results for i?(4)ij and -B(g)ij coincide 
with the ones from the last section. To see this we need 
the expressions for T(i)ij and T{e)ij- These follow from 
(3.9) after the transition to NW variables by (3.101 and 



PN expansion as 



'^6)u = E 



PaiSa(j)(k)^a + Paj Sa{i)(k)Sa 



' 2^^"'^"^^" ^ 2^^('')(^'^m^'' ^ 2Sa(k){J)n(l)5a + ■^^PaiSa{k)(3)<l^(2),k^a 
+ ^^PaoSa{k){i)<t)(2)M5a + ^^PakSa(k){i)<f)(2),j5a + ^^PakSa(k){j)(l>{2),tSa 



E A^' 

^ 4m3 



(5.32) 



(5.33) 



P'iPatSa{k){j}Sa + ^IPaj Sa{k){i)Sa + PaiPakPalSa{l)(j)Sa + Paj PakPal Sa{l){i)Sa 



For our source it holds 7(4);; — 2n'^^^*^°'^ and thus A ^T(i)ii — — g^0i(4). The quantities Bi^^^ij and B^^g-fij, 
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and thus the evolution equations of the gravitational 
field, can now be seen to coincide with the result of the 
last section. 



VI. ENERGY FLUX 

In this section we reproduce, within the ADMTT 
gauge, the IPN energy flux at the SO level, obtained in 
[1 , 12 within the harmonic gauge. This should be seen as 
a further check for the wave equation (with source terms 
depending on standard canonical variables) derived in the 
last section, which is most important for the calculation 
of the 3.5PN Hamiltonian. It also gives a check of the 
applied regularization techniques. Notice that the New- 
tonian flux in the SO and S1S2 cases vanishes. We will 
restrict to two objects here. 



A. Far zone expansion of the wave equation 

The retarded solution of the wave equation reads 

i3 //(^'j^ret) 



□ -J/(x,t)^-i- J 



(6.1) 



where / is some field and trot = t — c^^|x — x'|. We ex- 
plicitly show the speed of light c here in order to simplify 
the discussion. / shall not change much during time in- 
tervals significantly smaller than some time interval T, 
e.g., / describes a binary system with an orbital period 
T and the internal dynamics of the individual objects 
does not introduce a significantly smaller time scale rel- 
evant for /. This allows Taylor expansion in time of / 
in certain cases. The first case is the near zone defined 
by |x — x'l <^ cT where /(x'jiret), and thus D~J^f, can 
formally be expanded in c~^. The near zone expansion is 
important for the calculation of the Hamiltonian, as the 
metric at the position of the spinning objects is needed 
there. The second case is the far zone (or wave zone) 
defined by |x'| < cT < |x| = i? (if the support of / is 
centered around the origin). Again we can formally ex- 
pand in , but the quantity t^lt — ^ ~ ^ must be held 
constant. Both near and far zone expansion thus fit well 
into the PN scheme, as they can be seen as expansions 
in c^^. Useful formulas for the far zone expansion are 



|x-x'| = R-n-n' + 0{R-^ 

1 _ 1 
|x-x'| ^ R 

R 



R 



^0{R-^), 
fiK'Xl,)+OiR-^). 



(6.2) 
(6.3) 

(6.4) 



where = x'^ /R and L is a multi- index with / even. 
Notice that the spacial derivative of t^^t does not vanish. 



Equation (6.4 1 can be shown using 
1 



A-'sR-'e' 



R 



i-i 



3=0 



(ffi?v 



.5) 



and the Fourier transform of /(x',t^^j) with respect to 
time. 



(6.6) 



/(x',C)= J d^c./(x',c.)e-*-. 

This Fourier transform can also be used to show 
1 



dVe="'""''^'/(x',t::j+0(i?-2), 



(6.7) 

from which the far zone expansion in c ^ can most easily 
be obtained; one just needs to plug in the Taylor series 
of e'^ ^'^ ' to the required order. This is precisely the 
multipole expansion of the far zone field. 



Now we write the wave equation for ft/j , Eq. (|5.15|), 



in the form 



We can replace Ski by its STF part 



^kl^ — T^i^kl + Sik) — -5klSit 



here. The IPN far zone expansion of 
results in (from now on c is dropped again) 



(6.9) 



(6.10) 



'-kimn 



(6.11) 



+ OiR-^), 



with the multipole moments 



(6.12) 



Here M is a multi- index and the moments Ium are STF 
with respect to k and I. At higher orders it is better 
to work with multipole moments which are STF in all 
indices; see, e.g., 1501 HI] and references therein. The 
TT-projector Sj^''-' was replaced by Pijki, 



P, 



ijkl 



Pi(kPl)j ^ 2^ijPkl , 



Pij = Sij - n'n 



(6.13) 
(6.14) 



using (6.4). The integrations needed for the multipole 



moments also appear in the calculation of the Hamilto- 
nian. Details on the calculation and the applied regu- 
larization techniques will be given in |52| . we only show 
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the results here. Though the expressions for the multi- 
pole moments are quite long, after extracting total time 
derivatives as 



(6.15) 
(6.16) 
(6.17) 





— Qij 7 


lijk 


— Qijk 




— Qijkl 



they can be written in the compact form 



Qij — 



mizlzl ^—PikSuk)(i:)zi 

mi 



STFij 



+ (1^2), 



Q 



ijk 



iijkl 



2Piiz{zl - z\z{Pik - 2z'i5i(j)(fc) 

-(lo2), 
2 

TTii 

+ (1 o 2) 



(6.18) 



STF., (6.19) 



STF., (6.20) 



Here the subscript STFij means to take the STF part in 
i and j, and (1 O 2) denotes an exchange of the particle 
labels. Only terms needed for the IPN SO part of the 
flux are shown. As the source of the wave equation is ex- 
pressed in terms of variables with a standard canonical 
meaning, leading order Hamiltonians are used to calcu- 
late the time derivatives appearing here. 



B. IPN energy flux 

The energy flux C results from the formula 

lim (f dnhj^hj^ , (6.21) 



C = 



327r 

where is held constant in the limit R oo. Using 
(6.11 ) we can express this in terms of Iki, Ikim, and Ikimn 



as 



^IPN — ■r(-^ij)^ 



35e 



2 I o K-'-tjk) '^J-ijk-'-ikj 



'^{-^ikk} ~ ^-^ij i ikjk ^~ i ijkk 



3.22) 



A symbol (n) on top of a multipole denotes the n-th time 
derivative. Plugging in our expressions for the multipole 
moments we get for the SO part 



, , 12M \ 
27f^ - 37v2 — j 



P21 ( 18r2 - 19v2 - 



(6.23) 



+ (1^2), 



where Sa has components iSaj^jj, after going to the center- 
of-mass frame (see Appendix IdI) , where it holds 

M = mi + m2 , (6.24) 



P21 

v 

ri2 

r 

P 
L 



m2 
mi 



Pl2^ 



P2l 



Pl2 



Zi - Z2 , V 

||ri2||, ni2 

Pi =-P2, 

ri2 X p. 



(1+P12)2 
= Zi - Z2 

- ^ 

r 



(6.25) 

(6.26) 

(6.27) 

(6.28) 

(6.29) 
(6.30) 



Our result for Cso exactly coincides with the one in [1] 

Hj. 

Similarly, we could get the total angular momentum 
flux. This requires one to keep all O (R^^) terms. 



VII. CONCLUSIONS AND OUTLOOK 

In the present paper we extend the ADM canonical 
formalism to spinning objects up to and including 3.5PN 
and linear in the single spin variables. Further, general 
formulas for the interaction Hamiltonian and the wave 
equation for hj^^ are derived. This is the foundation 
for the calculation of the 3.5PN SO and S1S2 radiation- 
reaction Hamiltonians in |52j (which, respectively, are 
4PN and 4.5PN for maximally rotating black holes), as 
they result from the general interaction Hamiltonian by 
utilizing a near zone expansion of the wave equation .50] . 
The important difference to the formalism in |19j is the 
spin-dependent correction to the canonical field momen- 
tum. 

The conservative NNLO SO and S1S2 Hamiltonians 
(which are formally at 3PN or, respectively, at 3.5PN 
and 4PN for maximally rotating black holes) can also 
be obtained from the results of the present paper, once 
they have been rederived in arbitrary dimension. This is 
needed for dimensional regularization, which is the only 
one known to work consistently at 3PN in the ADM for- 
mahsm; see [S31[51]. The NNLO SO Hamiltonian is the 
last missing piece to complete the EOM for maximally 
rotating binary black holes up to and including 3.5PN. 

As the calculation of the mentioned NNLO Hamiltoni- 
ans is quite involved (comparable to the 3PN point-mass 
Hamiltonian), it is a good idea to thoroughly check the 
used formalism before starting such a calculation. The 
present paper provides several such checks. The wave 
equation for hj^"^ was compared with the Einstein equa- 
tions and applied to the leading order SO energy flux. 
Further, agreement with the action approach in [41 up 
to the considered order was shown. 

The method to construct a canonical formalism in the 
present paper is formulated in a quite general way and 
could be applied also in other situations. In particular, 
our method makes no use of a covariant generalization 
of flat-space expressions, in contrast to [19]. It is thus 
applicable to nonminimal couplings, which appear at the 
Si level. One can use the four-dimensional stress-energy 
tensor with Si quadrupole terms |551 156j to derive the 
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NLO Sf Hamiltonian with the method of the present pa- 
per. This Hamihonian was already obtained in [40l |55] , 
and agreement with the spin EOM from [57] was shown 
in [SS]. The LO conservative dynamics was obtained in 
[Tin]- For more LO resuhs at the level oriented toward 
application in GW astronomy, see, e.g., [59H63] . 
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Appendix A: PN orders and spin 

In the present paper, PN orders (orders in c~^) are 
counted in terms of the velocity of light c originally 
present in the Einstein equations. We call this formal 
counting. This has some computational advantages in 
our method, e.g., similarities to calculations for nonspin- 
ning objects are more manifest. Further, this formal way 
of counting best reflects the computational demands, e.g., 
the difficulty of the integrations and the regularization 
techniques that need to be applied. 

On the other hand, it makes sense to assume that the 
spin variables possess a numerical value of the order c~^, 
which holds for maximally rotating black holes. This way 
of counting best reflects the relevance of the spin correc- 
tions to the motion of rapidly rotating objects. Com- 
pared to the formal counting, it adds half a PN order 
for each spin variable appearing in a specific expression. 
For example, the NLO SO and S1S2 Hamiltonians are of 
the orders 2.5PN and 3PN for rapidly rotating objects, 
respectively, but are both of the order 2PN in the formal 
counting. 

To conclude, the formal counting overestimates the im- 
portance of the spin corrections, while the other way of 
counting overestimates the computational complexity. 



Appendix B: More on the Poincare algebra 

The action of an element of the three-dimensional Eu- 
clidean group, a subgroup of the Poincare group, on the 
coordinates of the three-dimensional hypersurfaces can 
be written as 



(Bl) 



with a* a constant infinitesimal vector describing a trans- 
lation and a rotation matrix Rij{uj). It holds R{uj) = e". 



where u)^^ = w-'* is a constant antisymmetric matrix de- 
scribing the axis and angle of the rotation. This is the 
standard representation of the Euclidean group on the 
coordinates. On a field, e.g., the metric, the standard 
representation of the Euclidean groups acts as 



7jj(x) Rik{uj)Rji{uj)"fki{R ^(w)x- a) 



(B2) 



Obviously, the ADMTT gauge conditions (2.7) mani 



festly respect the Euclidean group in its standard rep- 
resentation. Thus the global Euclidean group, as a part 
of the global Poincare group, is given by its standard 
representation in the ADMTT gauge. 

Now we restrict to infinitesimal transformations, i.e., 
a* and w'-' shall be small. Then it holds 



x'' ~> x'' + Sx'^ , 



(B3a) 
(B3b) 



where M, 



kl 



8^5i. The matrices M^- with compo- 



nents {Mij)'^^ = M^^ satisfy the commutation relations 

[Nhj , Mm] = 5^kMJl ~ 6jkMu + SuMkj - SjiM^ (B4) 

(B5) 



Thus the matrices Alij form a representation of the Lie- 
algebra so(3), namely, the vector representation. Further 
from (|B2|) we have 



7y ^ 1^3 - l^'^^M'ki 
which can be written as 



(B6) 



(B7) 



where C denotes the Lie-derivative. [This would of course 
be valid for any infinitesimal coordinate transformation. 



not only for (B3| 



As the Euclidean group is given by its standard repre- 
sentation in the ADMTT gauge, its generators in phase 
space are also given by its usual representations, Eqs. 
([Z25l), ([2^, and (fSlOl). Indeed, the transfor- 



mation rule for an arbitrary phase space function A, 

A^A+ Ilu'^A, J, J + a'{A, PJ , (B8) 
applied to our fundamental variables then reads 

z:^z: + iL.'='M2]x^+a% (B9a) 

Pa^ ^ Pa^ + i^'^'M^^P,, , (B9b) 
Sa(t)U) Sa(t){j) + \uj''''{Ml"i''Sa(m){j) + M^"S'a(i)(„) ) , 

(B9c) 
(B9d) 
(B9e) 



ijTT 



Thus the generators and Pi in its standard represen- 
tation give the transformation induced by ( B3 ) on the 
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fundamental variables, as expected 
it was used that, e.g., ^^^'^^ 



In (B9dl and (B9e 
6™ Cs.^hJ^ = Cs...hj;^ which 
again reflects the compatibility of the ADMTT gauge 
with the standard representation of the Euclidean group. 

Further, the canonical action given by Eq. (51) in [IT] 
(see also (4.33) in [TO]) is invariant under 



kl 



tj^[i](m) 



(BIO) 



and the transformations (B9). The corresponding 



conserved quantities can be obtained in the standard 
Nocther manner and result as (2.25) and (2.261 with 



(2.29) and (2.30) inserted, as expected 



A straightforward (3+l)-split of the Poincare algebra 
leads to 



{E,P^} 
{ Jfci, Jji} 

{G^p^} 

{G\G^} 
{G\E} 



0, 
0, 
0, 



= -J, 



I] , 



p 



(Bll) 
(B12) 
(B13) 
(B14) 
(B15) 
(B16) 
(B17) 
(B18) 
(B19) 



In consideration of ( B8 1 the first two equations reflect 



the translation invariance of Pj and E, while the third 
one requires to be a scalar under rotations. Simi- 
larly, the next equations state that Pk and G*^ trans- 
form as vectors under rotations, while Jki transforms as a 



bivector. Equation (B17) means that the center-of-mass 
j^k Qk jjTj j^g^g ^]^g expected transformation property 



under translations, {X^,Pj} 



Thus all except the 



last two equations are fulfilled by construction if Jij and 
Pi are given by its standard representation"^. For the 
calculation in 3D] the fulfillment of (B19) implied that 



(B18) also holds. However, a generalization of this fact 



is not known to the authors. 

It is instructive to give a physical interpretation of 
(B18) and (B19). Equation (B19) can be written as 



X' ^ {X\E} = ^ ^ const, 
E 



(B20) 



and states that the center-of-mass is moving with con- 
stant velocity. If we define a total spin of the system 
as 



-rtotal 



J j. 



X'P, +X^P, 



(B21) 



we get 



{X\X^} = 



Ctotal 

E-^ 



r ctotal Tj \ 



|5,total^^fc| 

r ctotal ctotal T 
I'-'ij ' "^kl S 



0, 

p ctotal 

E^ 

Skm 



p ctotal 
^j'^ik 



E^ 



PkPrr 



rnin ctotal 



Sim — 



E^ 
PlPrn 



E^ 



(B22) 
(B23) 

(B24) 
(B25) 



mn rftotal 



kn 



These are the Poisson brackets known for the center and 
spin associated with the SSC S^f^^^ = 0. Notice that 

^lf''\E}^0. (B26) 

One can go over to NW variables by 



ctotal 

ij 



total 



TD I ) ctotal 

P.PkS,^ , r-,r-,o,,^ ^^27) 



ototal 

"'^ ' M{E + M) ' M{E + M) 

^ total 

= X'- 



^total^. 

ij J 



M{E + M) ' 

where IvP — E^ — P^; see, e.g. 
(IMTl) into 



(B28) 
64l. This transforms 



s: 



total 



, (B29) 
and finally leads to the standard Poisson brackets 

{Z\P,} = %, (B30) 
{Z\Z^}=0, (B31) 



l^total^ ^fe| 
total ototal 



0, 



(B32) 



} 



M^J-Sl^r' + MirSl°:r'- (B33) 



It still holds Z^ ^ P,/E ^ const and S 



total 



0. Equa- 



tion (B18) was transformed into (B31). Thus (BIB) re- 



fiects the fact that by a straightforward (3-|-l)-split of 
the Poincare algebra one arrives at a center X* associ- 
ated with the (noncanonical) SSC 5'^°*''' 0. 

Notice that the structure of the total angular momen- 
tum in (B21 ) and (B29) is the same, but the variables in 



(B21) are not standard canonical. It is thus astonishing 



that the condition ( 2.30 ) uniquely fixes the canonical spin 
and position variables in this paper. This is due to the 
important additional requirement of having a constant 
Euclidean spin-length. In this section, however, even the 



individual components of 5*°'^^ and 5'*°*''' are constant 
due to the fact that the total system does not interact, 
e.g., with an external field. 



Appendix C: Lie-shift version of the transition to 
the NW position variable 



^ If is determined by an ansatz instead of the integral ( |2.24| 
then one should also check (B17l. 



For Szl^ = one can use a Lie-shift to redefine the 

position in ^matter j^j^j ^matter ^ j 

-^matter ^ -^matter -^matter (CI) 
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^matter ^ ^matter -^matter (^2) 

with the shift Sx^ on the a-th worldhne given by 



npa 



(C3) 



It holds Sx'^ — and Pai5x\ = 0. Although spatial 
derivatives of Pai are not defined (the linear momentum 
is only known on the worldline) , Pai is treated as a vector 
field for the Lie-shift. Therefore we must have 



= 0. As Pi is conserved 



const (cf. Sec. II B| this 



mass vector vanish, i.e., Pi = 

and = P,t + with 

is indeed a consistent set of constraints that can be im- 
posed on the phase space for all times t. We will restrict 
to two objects here. Then Pi — Q results in Pi ~ — P2 
at the considered order. The leading order terms of G' 
read 



G = miZi + TO2Z2 



Pi X Si P2 X S2 



2toi 2m2 



(Dl) 



5-^ 



^'a-iPak = -5x^a{Pai;k+Pak,i " Pai,k) = , (C4) 



which precisely cancels the spatial derivatives of Pai in- 
troduced by the Lie-shift. Thus pai^ as a vector field, 
must be parallel transported to the new worldline with- 
out rotation. The transition to the canonical momentum 
now has to read 



1 



Pai = Pai - nS^Kik - '^jk,i + -^Sajk 



p p 

am-'- aq 



nPa{ma - nPa) 



(C5) 



in order to satisfy Eq. (2.35[). For tt*-' = the momentum 



redefinition fro m [T^ is o btaine d. Further, this leads 
to the results (3.11) and (3.12) for the case ~ 0. 
However, this formulation is not so useful for variable 



redefinitions in an action approach. Notice that (C5) is 



missing a term when compared to (3.10c), however, here 
Pai is treated as a vector field and is not held constant 
for the redefinition of the position. 



From G = follows 



^ p X Si p X S2 

Zi — — ri2 , (D2) 

TOi 2miAf 2TO2M ' ^ ' 



with /i — vM, and similar for The leading order 
relation between canonical momentum and velocity reads 



Pu f 3to2 



mi .\2m] 



Si 



{i){3) 



and similar for z\. It follows 



25'. 



P ni2 X Si 



V = 



7712 
77ll 



2(oa) 



ni2 X S2 



''12 



(D3) 



777l\ 

7712/ 

(D4) 



The reduced phase space, where Zi, Z2, Pi, and P2 are 
replaced by ri2 and p, is then subject to the Poisson 
bracket 



Appendix D: Center-of-mass frame 



(D5) 



The center-of-mass frame is defined here by the condi- 
tion that the total linear momentum and the center-of- 



This can also be seen as a Dirac-bracket following from 
R- = G'' = 0. 
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